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Abstract 

i-£h . Decoupling of a heavy flavour in QCD is discussed in a pedagogical way. First 

ii | we consider a simpler case: decoupling of muons in QED. All calculations are done 

^H" up to 2 loops. 

^ 1 Introduction 

> ' 

QCD with all 6 flavours is rarely used. If characteristic momenta Pi <C Mq (where Mq is 
the mass of a heavy flavour Q), it is better to use a low-energy effective theory without Q. 
Its Lagrangian has the QCD form plus 1/Mq corrections. Operators in the full QCD are 
. expanded in 1/Mq via appropriate operators in the effective theory. 

! The pioneering paper [T] discussed decoupling effects in the MS scheme at two loops; 

however, it contained a calculational error. The correct two-loop result was obtained 
in [2] as a by-product of a three-loop calculation. A simple and efficient method to find 
decoupling effects was developed in [3]; the two- loop result [2] was confirmed, and new 
three-loop results were derived. After that, the erratum to [T] appeared, in which the 
authors fixed their error, and confirmed the results of [21 [3]. A few years ago, four- loop 
decoupling coefficients have been calculated [U [5jE 

In this lecture, we shall discuss decoupling effects in QCD at two-loop level. First we 
shall discuss the QED case which is very similar to QCD but simpler (Sect. |2]). We consider 
QED with electrons and muons; when characteristic momenta pi <C M (where M is the 
muon mass), the effective low-energy theory containing only electrons and photons can 
be used instead. After that we discuss QCD with light flavours qt and a heavy flavour Q 
(Sect. E]). 



* Lecture at 5-th Helmholtz international summer school Calculations for modern and future colliders, 
Dubna, July 23 - August 2, 2012. 

1 The result of [5] contains one master integral which was not known analytically, only numerically, with 
37-digits precision. An analytical expression for this integral has been published later [B]. 



1 



2 Decoupling in QED 



2.1 Full theory and effective low-energy theory 

Let's consider QED with nj = 2 lepton flavours — massless electron ip and heavy muon 

L = Vw'^oV'o + (ipo ~ M ) - -F^ V F^ - -L {d^f (2.1) 

(the covariant gauge is used; means bare quantities). When characteristic momenta 
Pi M (characteristic distances ^> 1/M), the low-energy effective theory containing only 
light fields can be used instead. Its Lagrangian contains all possible operators allowed by 
symmetries; those with dimensionalities > 4 are suppressed by powers of 1/M: 

L' = VM> - \F^Fr ~ (d.A'tf + C i°'i + ■■ ( 2 - 2 ) 

(quantities in the effective theory are denoted by primes). Muons only exist in loops of size 
~ 1/M producing local interactions of light fields. Coefficients in the effective Lagrangian 
are tuned to reproduce scattering amplitudes of the full theory expanded in Pi/M up to 
some order. We shall not discuss 1/M n corrections here. 

Operators of the full theory can be expressed via all operators of the effective theory 
which are allowed by symmetries. Contributions of higher- dimensionality operators are 
suppressed by powers of 1/M. Coefficients are tuned to reproduce on-shell matrix elements 
of the full-theory operators expanded in Pi/M up to some order. In particular, the light 
fields of the full QED can be written as 

i° 1 (2-3) 





We shall not discuss 1/M n corrections. 

Similarly, the parameters of the full-theory Lagrangian are related to those of the 
effective-theory one: 

eo=[C°]" V2 e , a =[a]- 1 a' . (2.4) 
We shall soon see why decoupling of the gauge fixing parameter a is determined by the 
same coefficient (a as that of the photon field. 

Now we shall recall some information about renormalization of QED which will be used 
later. For more details, see textbooks, e.g., [7]. 

The photon propagator has the structure 



MW\At = »WW» _J_ »V/\/V\/\^^^\/\/\/\/\» _J_ •WW\^^WW\^^WW\I _|_ . . . 

-iD, v {p) = - zD%(p) + (-i)Dl a (p)m^(p)(-i)Dl(p) ^ 

+ (-i)Dj a (p)zn^(p)(-z)D° 7 ( P )zn^(p)(-z)L>°>) + • • • 
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where the photon self energy ili^p) is the sum of all one particle irreducible diagrams 
(which cannot be cut into two disconnected pieces by cutting a single photon line). Due 
to the Ward identity, it is transverse: 

iWp) = (pV - p»Pv) n (p 2 ) , (2-6) 

and we obtain 

n f \ 1 ( Pi&A , VvVv , 7 x 

DM = ? [i - a ° w ' ( } 

The electron propagator has the structure 



iS(p) = iSoip) + iS (p)(-i)T,(p)iS (p) 

+ iSo( P )(-iM P )iS (p)(-iMp)iS (p) + 



(2.8) 



where the electron self energy — iE(p) is the sum of all one particle irreducible diagrams 
(which cannot be cut into two disconnected pieces by cutting a single electron line). For 
the massless electron, it has one Dirac structure 

E(p) = ^ v (p 2 ) , (2.9) 

due to chirality conservation. We obtain 

s M=w=km- (2 ' 10) 

We shall use two renormalization schemes, MS and on-shell. In the MS scheme, bare 
fields and parameters are related to renormalized ones as 



A = Z A /2 (a{fi)) A(fi) , ip = Z^ /2 (a{fi), a(fi)) ip{fi) , 
e o = Z l J 2 (a(ji)) e(n) , a = Z A (a(jj)) a(n) , 



(2.11) 



where all renormalization constants have the minimal structure 



*' = 1 + 7S + (7 + 7)(5) ,+ " < 2 ' 12 > 

(zeroth and positive powers of e are not allowed). We shall see in a moment why renor- 
malization of a is determined by the same constant Za as that of the photon field. In 
d = 4 — 2e dimensions, dimensionality of e is e; a should be dimensionless, therefore, we 
have to introduce renormalization scale fi to construct 

^ = ^^L-. (2.13) 
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Multiplying the bare photon propag ator flUD by Z A X converts D ± (p 2 ) to Dl(p 2 ; //) 
which is finite at s — > 0; at the same time it converts to a(/z) which is also finite. This is 
the reason why renormalization of the photon field A and the gauge parameter a is given 
by a single renormalization constant Za- 

The sum of one-particle-irreducible vertex diagrams (not including the external prop- 
agators) is the vertex function 



ie V"{p,p'). (2.14) 



It should be equal to T fl (p,p') = ZyT^(p,p'; fi) where Z? is a minimal renormalization 

1 /2 

constant and T^(p,p';fi) is finite at e — > 0. The scattering amplitude e T^Z^Z A = 
eT^Za 2 Z-pZ^Z^ 2 must be finite; this means that the product Z]j 2 Zy-Z^Z 1 ^ 2 is finite. 
But the only minimal (I2.12p renormalization constant finite at e — > is 1. Therefore, 
Za ZyZ^Z 1 / 2 = 1, and 

Z OL = \Z v Z i \~ 2 Z A X . (2.15) 

The Ward identity says that 




r*W)?M = s (2.16) 

Multiplying this equation by Z^ make its right-hand side finite; hence its left-hand side 
is finite, too, and Z? = ZT 1 . Thus the Ward identity makes situation in QED simple: 



Za 1 



In the on-shell renormalization scheme 

A = [Z° A s (e )} 1/2 A os , ^0 = [Z; s (e , a )] 1/2 ^ 



OS ) 



(2.17) 

e = [Z° s (e )] 1/2 e os , a = Z^(e ) a os 

(Z° s are not minimal). By definition, the renormalized propagators in this scheme tend to 
the free ones near their mass shells: 

DlV) DW) = 1 , S QS (P) 5b(p) = J (2-18) 

pr p 



at p — > 0. This means 



1 

1-Tl(0)' "W-u»«uy - !_ Sy ( ) • 



Z T( e o) = i ttTTvT . a ) = n ^ , n , • (2.19) 



When p and p' are on the mass shell, and the initial electron and the final one have 
physical polarizations, a single-photon scattering of an electron is described by 2 form 
factors; only one survives at q — » 0. The scattering amplitude in this limit is e os 7 M = 



4 



eoT^Z^ [Z^ s (eo)] 1//2 , where e os is, by definition, the on-shell electron charge. The vertex is 
T M = Zp s 7 M , and therefore 

ZZ=[Z?Z?]-*[ZZ\- 1 . (2.20) 
The Ward identity (12T6]) at q -»■ 

as- 1 (p) 



r^(p,p) 



near the mass shell (p 2 — >■ 0) gives 

r"(p,p) = 



<9p. 



L j* J 



(2.21) 



('2.221 



and we obtain Zf = [Z° s ] hence Z° s = [Z^\~ l . 

The MS renormalized fields and parameters in the full theory and the effective one are 
related by 

where the renormalized decoupling coefficients are 



(2.23) 



Z A (q(//)) 
Z A (a'(//)) 



^("(^))_ r 



They satisfy the renormalization group equations 

d logOiGu) 



(2.24) 



d log n 
d logCv-(^) 
d log fj, 

d logCa(yU) 



7a(«(m)) -7aO O)) , 
7^(a(/i), a(//)) - ^(a'(//), a'(/i)) 
2 - , 



(2.25) 



ti log [A 

where the anomalous dimensions and the /3 functions are defined as 

d log Z A {a(jj)) 



d log \x 
d log Zip(a(fi),a(fi)) 

d log fj, 
d \ogZ a (a(/i)) 



7a(«(//)) , 

7^(a(/i),a(/x)) 



(2.26) 



2P{a{n)) . 



d log /i 

Note that the MS charge a(/i) and the gauge parameter a(/x) satisfy the RG equations 
dloga(fj) . . dloga(/i) 



d log /i 



-2e - 2/3(a(//)) 



d log /i 



-7a (a(//)) • 



(2.27) 
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A os = A' os + 0[^-) . (2.29) 



2.2 Photon field and electron charge 

The propagators of both A os and A' os are equal to the free propagator at p 2 — > 0: 

D°l{p 2 ) = DT(P 2 ) [l + 0{p 2 )] , (2.28 

and therefore 

A^ = A' a- CI I 

,M 2 

Hence the bare decoupling coefficient is 

« = W (230) 

We have 

But n'(0) = because all loop integrals are scale-free. Hence Z'^ = 1: 

C = 1 - n(0) . (2.32) 
At one loop, the photon self energy f)2.6p is 

k + p 




. -„ / I /' </,,./ — /',,/>,, ! UK'"' 

jU \ / V 



k 

Contraction in /x and z/ we get 

4ie^ f d d k N 



n(V) =- 
{P ] {d-l)p 2 J (2tc)*D 1 D 2 

D 1 = Ml - (k + p) 2 , D 2 = Ml - k 2 , 
N = ^Tt + i> + M ) 7 ^ + M ) . 

Now we expand the integrand in p up to p 2 ; the problem reduces to trivial 1-loop vacuum 
integrals: 

n<o) = ■ (2 - 33) 

The bare decoupling coefficient of the photon field with the 1-loop accuracy is 



Re-expressing it via renormalized quantities in the full theory, we obtain 

= 1 + l^Z Q (a(ri)Z-*%a(ri) L^£j 2£ r(l + e)e* + ■■■ (2.35) 

where the MS renormalized muon mass is defined by 

M = Z m (a(ji)) M(ji) . (2.36) 

The MS renormalization constants of the photon fields in the full theory and the effective 
one with the 1-loop accuracy are 

Z^(a) = l-^n f £- + ... (2.37) 

(see, e.g., the textbook [7]). There are 2 lepton flavours, electron and muon, in the full 
theory (n/ = 2) and only 1 (electron) in the low-energy effective theory (n/ = 1). We 
may neglect the difference between a'(/j,) and a(fi) in corrections. Combining these pieces 
together, we arrive at the renormalized decoupling coefficient of the photon field with the 
1-loop accuracy 

( A (V) = 1 + ~L?P- + ... (2.38) 

3 47T 



where 



Note that L depends on /i in a complicated way because the MS renormalized muon mass 
M(/i) also depends on \x. It is convenient to do decoupling at /i = M where M is defined 
as the solution of the equation 

M(M) M; (2.40) 

then L = 0. The decoupling coefficient Ca(aO for other /x can be obtained by solving the 
RG equation (12.251) with this initial condition. 

The scattering amplitude of an on-shell electron with a physical polarization in an 
electromagnetic field at q — > should be the same in both theories: 

eo r^; s [z° A f /2 = e r'% os [zf] 1 ' 2 . (2.41) 

The vertex functions in the two theories are 

r M = z^ s Y , = z£ s Y , (2.42) 

and hence 

[ 70s 70s] 2 70s y/ os 

£ = 2 = ^- • (2.43) 

[Z^Z'^] 2 Zf Z£ 
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As we already discussed, Z^Z^ S = 1 due to the Ward identity; the same is true in the 

, ,.fl',„,t;,-„ (i„„„v 7'"" — \ 



effective theory. Moreover, in the effective theory Z^ s = Z'? s = Z'£ s = 1 because all loop 



integrals are scale- free. Therefore 

C-iClV 1 . (2-44) 

In other words, the on-shell electron charge (measured at large distances in smooth 
macroscopic fields, i.e., at q — >■ 0, in Millikan-like experiments) must be the same in both 
theories 

a os = «o S • (2.45) 

Thus we immediately obtain (I2.43p . 



Figure 1: The Green function of ipo, i>o, A)- 

We can look at the charge decoupling from a slightly different point of view. Let's 
consider the Green function of ^q, ipo, A . It is the vertex with the full propagators 
attached (Fig. [T]). It differs from the similar Green function in the effective theory by the 
bare decoupling constants of the three fields: 

e TSSD = [(J] - 1 [C] - 1/2 e' T'S'S'D' . (2.46) 

On the other hand, S = [C°] ~ l S', D = [Q' 1 D', and hence 

eo^ = Cl[Cl] ll2 e' T^. (2.47) 

Writing down r M = [Q]" 1 where Cr = Z' r os /Z° s , we obtain 

<2=[^]" a [CXI" 1 , (2-48) 

and this is equivalent to f)2.43p . 

The MS renormalized charge decoupling is 

= CVKM , (2-49) 

where 

<M = ^rr^- (2 - 50) 

Z' Q (a'(fi)) 

Taking into account = [Z^] , we obtain 

(M = Q 1 ^)- (2.51) 
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Figure 2: Two- loop diagrams for the photon self energy. 



With the 1-loop accuracy Ca(aO is given by (|2.38p . 

The photon self energy with the 2-loop accuracy (Fig. [2]) is 



n(o) = - 



4 e\M^ 



Tie) 



3 (iir) d / 2 
2 (d - 4)(5d 2 - 33d + 34) ( e\M^ 



-2e 



3 d(d - 5) 

It reduces to the integrals (Fig. [3]) 
d d h d d k 2 



\ (4vr) d / 2 



Tie) 



D n 1D n 2D n 3 



- 7r d M 2 ^- ni - n2 - n ^V(n 1 ,n 2 , n 3 ) 



3 

2 ;„2 



+ 



D 1 = M - k{ 



Do 



M — k 



2 i 



k,) 2 



(2.52) 



^2.53) 



which can be expressed via T functions [8] 

T (f - n 3 ) r (m + n 3 - § ) r (n 2 + n 3 - I) r(ni + n 2 + n 3 - d) 



V r (n 1 ,n 2 ,7i 3 ) 



T (|) r(ni)r(n 2 )r(n 1 + n 2 + 2n 3 - d) 



(2.54) 




Figure 3: Two-loop vacuum integrals. 
The bare decoupling coefficient of the photon field with the 2-loop accuracy is 



Ae 2 M^ 



Tie) 



3 (4vr) rf / 2 
2 (d - 4)(5d 2 - 33d + 34) /e 2 M, 



d(d - 5) 



-2s 



( 4vr )d/2 



■rfe: 



(2.55) 
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Re-expressing in via renormalized quantities in the full theory, we obtain 

= 1 + ^^Z a (a(/x))Z-*(aM) 



3 J \ Aire 

We need 1-loop corrections to Z Q and Z m in the full theory: 



(2.56) 



Z„ = Z^ 1 = 1 + 2 ■ --^ + ■ • • Z m = 1 - 3-^ + • • ■ (2.57) 
3 Ane Aire 

The MS renormalization constants of the photon fields in the full theory and the effective 
one with the 2-loop accuracy are 

zV(a) = 1 - -n f — - 2en f ( —Y + ■■■ (2.58) 

(see, e.g., the textbook |7j). Substituting everything into Ca(aO (12.241) . we arrive at 



a 1/ x a / x 4 r q(u) / AT 13 \ fa(a) 

cm = a(/i) = i + + (^- 4L + 



An 



2 



+ ■•• (2.59) 



where L is defined in ( 12.391) . 

In particular, at \i = M (f2T40D (L = 0) 

C«(M) = l-f(^)V- (2.60, 
If we use /i = M os instead, where the on-shell muon mass is related to the MS one by 

^)-l_ 6 ( r iog^- + ?U + ... (2.61) 



M os V. Mos 37 4vr 

then L = 8-p- and 

4-7T 

Ca(M os ) = 1 - 15 (j^X + • " • (2.62) 

In general, for any \i = M(l + 0(a)) we have Ca(A*) = 1 + 0{a 2 ) (with different coefficients 
of a 2 ); however, for, say, fi = 2M or M/2 the 1-loop term appears. It is most convenient 
to use some /j,q = M(l + 0(a)) as an initial condition, and obtain C«(aO f° r other \x from 
the RG equation fT2T23|) . 
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^ os = + O [ ^= ) . (2.64) 



2.3 Electron field and mass 

The propagators of both i/j os and ip' os are equal to the free propagator at p 2 — > 0: 

tS os (p) = ^S' os (p) [1 + 0(p 2 )] , (2.63) 

and therefore 

nh- = nh' 4- D ( 

Hence the bare decoupling coefficient is 

>7 Jos (J \ 

«=w (2 - 65) 

where 

Only diagrams with muon loops contribute to £y(0); £y(0) = because all loop integrals 
are scale- free. Hence Z'^ s = 1: 

d = l-E v (0). (2.67) 




Figure 4: Electron self energy. 
The first diagram contributing to £y(0) appears at 2 loops (Fig. H|): 

-^SvV) = y Al-i^Yi ieoY (j^j i(k 2 g^ - k^k v )Tl{k 2 ) , 

where U(k 2 ) is the muon-loop contribution to the photon self energy. Expanding the 
right-hand side up to linear terms in p, we obtain 

2 (d-l)(d-4) f d d k U(k 2 ) 



E V (0) = g i y . (2.68) 



The calculation reduces to the integrals (12.54)) . The result is 

_ 2(d-l)(d-4)(d-6) / egM - 2E , , 
M0) " d(d - 2)(d - 5)(d - 7) lWF^ ( } ' 
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+■ ■ 



The bare decoupling coefficient is 

2(d-l)(d-4)(d-6) ( e 2 M^ V J,_ 5 P , V°V 

d(d-2)(d-5)(d-7) U^/ 2 l 7 V 6 £ + "7U7r £ J 

(2.70) 

To find the renormalized decoupling coefficient (^(fi) ( 12. 24ft we need the MS renormaliza- 
tion constants Z$ . They are determined by the anomalous dimensions 

7?(«. «) = 2a£ - (4n, + 3) {^f + - ■ ■ (2.71) 

(see, e.g., the textbook [7J). Taking into account ac'(M) = a(M) [1 + 0(a 2 )}, a'(M) = 
a(M) [1 + C(a 2 )}, we see that 

Z^(a(M),a(M)) / a \ 



Z;(a'(M),a'(M)) V4vT£ 

and hence 



J 



^ = 1 + 5^'] + -. (2.72) 



6 V 4tt 

(Cv'Ca 4 ) f° r an y H ~ M has the same form, because there is no 0(a s ) term). To find Ci/>ip) 
for other values of /i, the RG equation (I2.25P can be used. 

Until now we considered QED with massless electrons and heavy muons. Now let's take 
the small electron mass m into account. The electron self energy has two Dirac structures 

Z(p)=^ v (p 2 )+m i;s(p 2 ), (2.73) 

and the propagator is 

S{P) = ji-mo-ECp) = 1 - MP 2 ) ~, 1 + S g (p 2 ) • (2 ' 74) 

1 - Z v {p 2 ) 



Near the mass shell 

i i wr 1 



1-E v (p 2 ) 1 + S s (p 2 ) "1-SV(p 2 ), 1 + Z' s (p 2 ) , 



(2.75) 



We shall work in the linear approximation in m. Comparing the overall factors, we repro- 
duce ( 12.65P ; comparing the denominators, we obtain 

1 + Ss(0) 1 + S'g(O) , 

I-Mo) mo = T^(o) mo (2 - 76) 
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(we may set itlq = in Ey j( sf(0)). The bare masses in the two theories are related by 

m =[C° m ]- 1 m' ; (2.77) 

we obtain 

,o = ™ -i 1 + S5(0) = 1 + Ss(0) r2 ^ 

w L ^ J 1 + S' S (0) 1-Ev(0) 1 J 

(because S^(0) = E' s (0) = 0). 

In other words, the on-shell electron mass is the same in both theories: 

m os = m' os , (2.79) 

because it can be directly measured in experiment. In view of m = Z™m OB and a similar 
relation in the effective theory, this leads to 

Z' os (p' \ 

C ° = ^Lpl . (2.80) 
Z«(eo) 

The on-shell electron mass renormalization constant in the full theory Z£f depends on two 
masses, m os and M os ; if we neglect corrections suppressed by powers of m^/M^, this (more 
physical) definition of Cm coincides with our previous definition based on expansion in m 
up to linear terms. 

The first diagram contributing to £s(0) appears at 2 loops (Fig. HJ). In the linear 
approximation, we retain mo in the numerator of the electron propagator, and set m = 
in other places. The result is 

Substituting it and (12.691) into (I2.78p . we obtain 

ao _ i _ 8(d-l)(d-6) ( elM^ \ 
W d(d - 2) (d - 5) (d - 7) V (47r) d / 2 1 7 



2 

+ 



5 89 2 \/a\ 2 



(2.82) 



To find the renormalized decoupling coefficient 



= #7^C (2-83) 

we need the MS renormalization constants Zm ■ They are determined by the anomalous 
dimensions 

7<>) = 6^+(3-H° ImV-. (2.84) 



4tt V 3 / / V4tt 
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(see, e.g., the textbook [7]); we obtain 



ZmWM)) = x _ ( 2 _ 5 \ (_a_\ 
Z'Ja'(M)) V 3 J We) 



2 

+ 



and hence 

, — . 89 fa(M)\ 2 

UM) = l + -(-^J +■•• (2.85) 

To find Cmil 1 ) f° r other values of /i, the RG equation 

C? log Cm (A*) 



<ilog /i 

can be used. 



+ 7m (a(/x))-7™(« / (/^))=0 (2.86) 



3 Decoupling in QCD 

3.1 Full theory and effective low-energy theory 

Now we shall consider QCD with n\ massless flavours and a single heavy flavour Q: 

m 11 
L = V q w ilpoqoi + Qa - M ) Qo - t^G ^ - — {d,A a ^f + (fygg) (Dffcg) . (3.1) 

When characteristic momenta <C M, the low-energy effective theory containing only 
light fields can be used instead: 

ni 1 1 / i \ 

L' = £ ?o^o4 - i ^G' ^ - — (9^D 2 + ( (itftf) + O (-) . (3.2) 

The fields and the parameters of the full theory are related to those of the effective theory: 

Ao = [C° A ] ^ A' , go = [Q " 1/2 g' , co = [C°] ' 1/2 c , 3 
9o=[Cr 1/2 9 f , a =[C]~ 1 a' . 
The gluon self energy up to 2 loops (Fig. E]) is 

n (°) - "3^74^7^ ( ) 



1 



d(d-5) 



2 /j ,w CJ 2 ooj d 5 -20d 4 + 145d 3 -458d 2 + 588d- 232 ' 
-(d - 4)(5* - 33rf + 34)<7, (d-2)(d-7) ^. 

-2 £ \ 2 



(3.4) 
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Figure 5: Gluon self energy up to 2 loops. 



(the Cp term trivially follows from QED; the calculation reduces to the vacuum inte- 
grals f)2.54p ). The bare gluon decoupling coefficient is = 1 — 11(0). 

The quark self energy up to 2 loops (Fig. HJ) can be trivially obtained from the QCD 
results 0231]), iEEU): 



= 2(d-l)(d-4)(d-6) c ^ {<,IV ( , 



s s (o) 



d(d - 2)(d - 5)(d - 7) 

2(d- l)(d- 6) 
~(d-2)(d-5)(d-7) 



(4vr; 



<f/2 



■r(e) + 



C,T, ( ^Cr(e) I i 



( 47r )d/2 



(3.5) 



Therefore, the bare decoupling coefficients for the light-quark fields (I2.67P and masses (I2.78P 
are 

2 



d(d-2)(d-5)(d-7) 

Q 8(rf-l)(d-6) 
W d(d-2)(d-5)(d-7) 



CfTf 



( 4vr )^/2 

^o 2 M - 2£ 
(47r) rf / 2 



(3.6) 



r(e) + 



The ghost propagator is 



G(p) 



1 



p 2 - X C Q» 2 ) ' 

therefore, the on-shell renormalization constant of the ghost field is 

1 



dYj c , , 
1-^(0) 



(3.7) 



(3i 



dp 2 

The first contribution appears at 2 loops (Fig. [6]): 



dp 2 



(0) 



2(d- l)(d- 6) 
'd(d-2)(d-5)(d- 7) 



CaTf 



(4vr)<V 2 



He) + 



(3.9) 
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Figure 6: Ghost self energy. 



3.2 Decoupling of a s 

All elementary vertices of QCD are determined by a single coupling constant go- Therefore, 
any vertex (quark-gluon, 3-guon, ghost-gluon, 4-gluon) can be used to obtain the decou- 
pling relation for g . In the full theory, we expand these vertex functions in their external 
momenta up to the first non-vanishing term, and compare them to the corresponding ele- 
mentary vertices in the effective theory (there are no loop corrections in it, because there is 
no scale) . This first non- vanishing term of the expansion ought to have the structure of the 
elementary QCD vertex, otherwise the Lagrangian of the effective theory would not have 
the QCD form (higher terms of expansions in external momenta lead to higher-dimensional 
operators in the effective low-energy Lagrangian, suppressed by powers of 1/M 2 ). 

The quark-gluon vertex function at 0-th order in its external momenta obviously has 
a single Dirac and colour structure / y fl t a . 

The 3-gluon vertex function should be expanded up to linear terms in its external mo- 
menta. The Bose symmetry allows two structures: / aia2(l3 (^ MlM2 (/c 1 — /c 2 )' 13 + cycle) and 
d a ia2 a s (^iw/jW + cycle ) The Slavnov-Taylor identity <T {d* A^x) , d u A u (y), d x A x (z)}> = 
leads to ^l\%%k^ k%* k% 3 = 0, thus excluding the second possibility. 

The ghost-gluon vertex has a single structure: p^f abc y where p M is the outgoing ghost 
momentum. We shall prove this statement in a moment. 

As a result, each vertex function of the full theory, expanded in its external momenta 
up to the first non-vanishing term, is equal to the corresponding elementary vertex times 
a scalar quantity IY Any one of these vertices can be used to find the bare decoupling 
coefficient (®: 



We shall use the ghost-gluon vertex here, because the calculation is simplest in this case. 

Now we shall discuss the ghost-gluon vertex function in the full theory expanded in its 
external momenta up to linear terms. Let's consider the right-most vertex on the ghost 
line: 



The tensor A^ v may be calculated at zero external momenta, hence A^ v = Ag^ v . Therefore 
all loop diagrams have the Lorentz structure of the tree vertex, as expected. 



(3.10) 
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Now let's consider the left-most vertex: 




It gives k , thus singling out the longitudinal part of the gluon propagator. Therefore, all 
loop corrections vanish in Landau gauge. Furthermore, diagrams with self-energy insertions 
into the left-most gluon propagator vanish in any covariant gauge: 



In the diagrams including a quark triangle, the contraction of k transforms the gluon 
propagator to a spin propagator and a factor k p which contracts the quark-gluon vertex. 
After decomposing ft into a difference of the involved fermion denominators one obtains in 
graphical form 



The diagrams with a massless triangle vanish. The non-vanishing diagrams contain the 
same Feynman integral, but differ by the order of the colour matrices along the quark line, 
thus leading to a commutator of two colour matrices. 

The remaining diagram contains a three-gluon vertex with a self energy inserted in the 
right-most gluon propagator. The contraction of k x with the three-gluon vertex cancels 
the gluon propagator to the right of the three-gluon vertex: 
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The colour structure of the three-gluon vertex is identical to the commutator above, how- 
ever with opposite sign. Therefore, after summing all contributions the complete 2-loop 
result is zero. 

Thus the bare decoupling coefficient £° with the 2-loop accuracy can be obtained (13.1 Oft 
from the on-shell field renormalization constants Z A S (13.4ft and Z° s H 3 . 9 1) : 



[C] 



+ 



d-4 
d(d - 5) 



V-3M + 34)g,- J - 1 ^ + f 



( g 2 M 

V (4vr) rf / 



-2: \2 

r(s) +■■ 



(3-11) 



It is gauge invariant, as expected. At 2 loops, Z A S and Z° s don't depend on ao; this means 
that Tacc does not depend on ao at this order, and its 2-loop term trivially vanishes in 
Landau gauge a = 0. Therefore non- vanishing 2-loop diagrams for Tacc must cancel, as 
discussed above. At 3 loops this is no longer so. 

We have the expression for g 2 = CaG?o)<?o y i a the bare quantities of the full theory. We 
re-express it via the renormalized quantities using 



( 4vr )<V: 



-V(e) 



'lire 



(3.12) 



where 



Z a (a) = 1 - A 



a 
4ire 



a 
Aire 



+ 



and (12T36|) . 

Inverting the series 



(3.13) 



9t 



( 4vr )d/2 



Tie) 



(3.14) 



we obtain 

47T£ 



/2 /-2e 
^0 /f c - 7 £ 



( 4vr )rf/2 



1 + ^f^ e_7£+ f /3 ° 2 + ^ 



/2 /-2e 
#0 /f - 7 e 



( 4vr )d/2 



(3.15) 



Finally, we substitute the expression for g 2 via a s (/x) obtained above, and arrive at a' s (/i') 
expressed via via a s (/i). 

The final result for // = is a' s (/i) = ^Q,(/x)a s (/x), where the renormalized decoupling 
coefficient is 



CM 
+ 



-_ i _ i LTF ^M 

3 4tt 
rpL 2 + 4 ( - -(', 1 L 



13 32 



Tt, 



47T 



(3.16) 



+ ■ 
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(L is defined by (12.391) ). It is convenient to use fi = M defined by (12.401) : 



(3-17) 



(here the Cp term trivially follows from QED (I2.60p ). For other values of /i the RG 
equation ( 12. 25ft can be used. For example, Fig. [7] (produced using RunDec [9]) shows 
af\/j,) and af 1 ^) near fi = Mb. 

0.225 



0.22 



0.215 - 




0.21 
M° s - 0.5GeV 



M° s + 0.5GeV 



Figure 7: Crossing the 6-quark threshold: as (/J,) and af\lj)- 

The decoupling relation for the light-quark masses m'(M) = ( m (M)m(M) can be 
trivially obtained from the QED result (12.851) : 

Um = 1 -*1 CF T F (^)\... (3,8) 
For other values of /i the RG equation (12.861) can be used. 



„(4), 



3.3 QCD fields 

Decoupling of the gluon field and the gauge parameter are given by the same quantity 

a 'o = (1(90,0-0, M )a . (3.19) 

At the first step we express the bare quantities in the right-hand side via the renormalized 
ones using (13.121) . (12. 36 p . and 



a = Z A (a s (/j), a(/i)) a(/i) 



(3.20) 
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and thus obtain an expression for a' via the renormalized parameters of the full theory 
a(fj,), M(/i). Next we find a'(/jf) in terms of a' by solving the equation 

a' = Z' A (a f M,aW)aW (3.21) 

iteratively. Substituting the expression for a' obtained earlier, we arrive at an expression 
for a'(fi') via a(fi), a s (/i), M(/i). 

It is convenient to use // = \i = M: a'(M) = ( A (M) a(M), 

U{M) = l+ l ^{AC F -C A )T F {^^ +•■■ (3.22) 

(the Cp term trivially follows from QED (I2.59P ). Coefficients of this expansion are gauge- 
dependent starting from o? s . For other values of fi, // the RG equation (12.251) can be 
used. 

For the quark field decoupling we can trivially obtain from the QED result (I2.72p 

(qi M ) = 1 + lcrT F (?4§l) 2 + - (3.23) 



6 V 47r 

Coefficients of this expansion are gauge-dependent starting from 



For the ghost field we find 



89 fa s (M) x 1 



( c (M) = l--C A T F l -±-L) +.-. (3.24) 



4 Conclusion 

Decoupling effects in QED are tiny; they are discussed in Sect. [2] for pedagogical reasons: 
QED calculations are similar to QCD but simpler. 

In QCD decoupling effects are essential for any quantity if we want to study it in a wide 
range of \i. The most important application is a s (ji): if we want, e.g., to relate a s {m T ) to 
a s (mz), we need to take into account both RG running (in each interval it is controlled by 
the corresponding (3 function) and decoupling. The same is true for, say, m s (fi). Parton 
distribution functions are also extracted from experiment at very different values of /z, and 
in addition to RG running (given by the evolution equations) we need to take into account 
also decoupling effects. These effects are more difficult to calculate than, e.g., for a s (/x); 
however, the methods used are similar. 

QCD without one or more heaviest flavours is one of the simplest examples of an 
effective low-energy theory. There are other similar examples, e. g., the effective Lagrangian 
of the Higgs-gluon interaction. In the Standard Model, it is produced by the t-quark loop; 
however, characteristic momenta are <C M t , and this loop can be replaced by a local 
interaction. 

I am grateful to K. G. Chetyrkin, M. Hoschele, J. Hoff, and M.Steinhauser for discus- 
sions of decoupling, and to D. I. Kazakov for inviting me to give a talk at Calc-2012. The 
work was supported by RFBR (grant 12-02-00 106-a) and by Russian Ministry of Education 
and Science. 
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